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Tee ROUT LON 


i) PROBLEM 

A problem nas occasionally been experienced in concen- 
trating laser energy from a laser designator onto a target 
and sufficiently illuminating that target so a weapon 
sensitive to that illumination will home in on the target. 
Two circumstances have possibly caused this to happen; the 
energy is either too diffused at the target to illuminate it 
properly or the beam partially or completely misses the 
Target. The aim of this study is to develop a computer 
program to aid in determining which part of these effects is 
atmospheric and which is the instability of the laser and 


Peaolliization system in the illuminating aircraft. 


Bee METHOD OF APPROACH 

The analysis of the problem requires two general types 
of systems. One system evaluates the performance of the 
designator using some techniques from Scott [Ref. 1]. A 
Bauteon vVidicon views the successive pulsed deSignator spots 
on the target from a location near the target. A line 
Spread function of the spot on the target and the standard 
deviation of the wander of the laser beam iS produced. The 
other system provides measurement of eae ior -avcmoespneric 


turbulence along an optical path adjoining the designator 





optical path. The laser on the target is viewed by a 
silicon vidicon located near the designator's optics. fhe 


- obtained is inherently properly path weighted 


value of C 
to express the atmospheric effects on the deSignator spot. 
Techniques from Fourier optics theory use the measured value 
of eae to predict the intensity pattern of the laser beam on 
mmemtistant target. The results from the two systems can then 
be analytically compared to determine the causes of spot 
wander and broadening. 

The sponsor-supplied hardware includes an HP-9825 
Calculator with 23K bytes of internal memory for equipment 
@emeuero! and data processing; an HP-9885 dise memory used for 
additional storage space; and a Quantex DS-30 Digital Video 
Analyzer used for digitizing the analog video input from an 
Eigen video disc; a Panasonic NV-1240 video tape recorder 
memeinacial recording of video from the target on a silicon 
feereon Camera. After analysis, output is in the form of 
merecs produced on an HP-9862 plotter. 

The sequence of analysis is currently as follows. 
Meaulation transfer functions (MTF) of the laser output and 
optics are measured or calculated and stored on dise. A 
video recording of the target being illuminated by the 
designator is made for a number of video frames, the target 


momen is also recorded, for later subtraction from the total 


picture to produce an image of the laser spot alone. The 





video recording is then played through the DS-30, controlled 
By the HP-9825, to digitize the information for use by the 
HP-9825. The HP-9825 takes the digitized video of the 
target, produces a line spread function (LSF) of the image, 
Subtracts the background, and averages a number of frames to 
produce a short term measurement of the laser energy on the 
target. The amount that each laser spot wanders from the 
others is also stored. In the second phase, the HP-9825 
uses the MTIF's of the laser, optics, and atmospheric 
measurements to produce a total system predicted MTF, LSF, 
and wander. These results of the target spot measurement 
and the baseline studies can be compared for possible 
Borrelation. 

Mieeelaset pulse 1S timed so that it occurs on the fly- 
back of the TV signal. This will cause the laser spot to 
appear on one field of aframe. For accurate results the 
BeeoO Needs to digitize a single field for each pulse of the 
laser. The videotape unit currently used will not allow 
mocecapablility. For this reason the tape recording video 
Meecnansferred to an Eigen video disc, which has the 
Sepablility of displaying each frame or field (1/2 video 
frame interlace) individually and of selecting any one frame 


Steer leld for viewing. 





II. BACKGROUND AND PROGRAM DESCRIPTION 


mee PEEORETICAL DESCRIPTION 

In order to combine the individual optical components of 
the laser to the target system, some elements of Fourier 
theory must first be remembered. The Frourier transform 
represents the one dimensional position variable g(y) 
expressed in the spatial frequency domain U(v). The inverse 
transform repeats the operation in the opposite direction. 
Their forms are commonly represented as follows: 


te] 


U(v)= , PA ex pce Teivy)dy= #igty) 
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gly = { Uy ex ole nr ivy) dv = F-'tuv)] 


imemprogram calculates the value of the integral at a preset 
number of points and yields a discrete Fourier transform. A 
Beeorem from Fourier tneory that allows us to compute the 
line spread function as presented in Fried [Ref. 2] for one 
of the elements in the system given the others is the con- 


volution theorem stated as follows: 


ine maneaey ice. G CY) 
anduasee (ncy) l= H(v) 
then &# [g*h] = GH 
(* - convoluted with) 





The first element in the system we are analyzing is, of 
Sourse, tne iaser itself. Ihe radial intensity output of 
the laser is generally Gaussian of the form Men exp (or -/ea-) 
as shown in Figure 1. This can be calculated for each laser 
Pereapproximacte results but was usually measured directly in 
mais Study. 

The above radial distribution of intensity must then be 
integrated over x to produce a one dimensional line spread 
mametLon for the source. A Fourler transform is then used 
Smeetis LSF to yield the MTF of the laser. The MTF is a 
measure of the spatial frequency response of a system 
compared to the input. A "perfect" system response would be 
Mmeorout to its limit (large spatial frequencies) and then 
Smopping to Zero. 

DG@emopeles associated with the laser is basically 
meerimaccion limited. The “Airy function" is used to cal- 
Pueeacve the diffraction point spread function for the laser 
meeles. Ihis point Spread function is then converted to a 
mene spread funetion, Fourier transformed to get the MTF and 
Chen combined by the convolution theorem into a laser system 
Bourlier transform. 

Mm@emnext step 1S to calculate the effects of atmospheric 
turbulence on the laser beam. The value of Ca Pine Vundex 
Sierefraction turbulence structure constant," as expressed 
Pyaeune relationship from Tatarski [{[Ref. 3] and Ochs et. al. 
(Ref. 4]: 
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where La es es is absolute temperature at two points 
Beperaced by a distance r and p iS pressure in millibars; is 
used to predict the MTF of the atmosphere over the path 
between the laser and the target. Measurements made along 
the path are inputs to the Ce equation which predicts the 
turbulence effects on the MTF and wander of the beam. In 
the field, however, plans are to measure on by optical 
methods decribed in Crittenden [Ref. 5]. 

The theory then, according to Fried [Ref. 2], shows that 
the total effect of the atmosphere and the laser can be 


aod bDY¥e tne product of the two transform functions: 
Cerourdle(ve-=G opvics(v) % G atm(v) 


The inverse transform then will yield the line spread 
function of the total system, the LSF as predicted on the 
target. 

Crittenden et. al, [Ref. 5] describe the process of 
converting the one dimensional LSF to a circular symmetric 
PSF using the Abel transform, from Griem [Ref. 6]. The 
power inside a circle of radius R can then be obtained by 
Maeeeratving the PSF out to ®. For the purposes of compari- 
son in this study, however, the resulting system LSF is used 


to compare the prediction with measured values. 
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B. HARDWARE USE AND LIMITATIONS 

(Meme =9325 Calculator used for this program had an 
internal capacity of 23K bytes. This was quite adequate for 
program storage but because of the great volume of digitized 
pcieo Gata that was involved, the inclusion of an HP-9885 
dise drive system was necessary. The Quantex DS-30 can 
Score and rapidly transfer a maximum of 512 lines of video 
with 512 pixels (picture elements) per line, easily 
exceeding the HP-9825's memory capability. The large 
overflow of data was reduced by storing each image into a 
quarter of the memory, thus limiting the amount of data that 
needed to be transferred. This also made disc storage more 
reasonable and sped up program running time. 

The DS-30 was also capable of taking the difference on a 
pixel by pixel basis between a reference image and the input 
and storing that in memory. Differencing was necessary in 
order to remove the background from the recorded image for 
processing the laser spot alone. Normally, it is necessary 
@Merwdivide the image by the background to offset the effect 
memone nOn—uniform reflectivity of the target. This was not 
possible with the DS-30 so these steps need to be performed 
Dy the computer. The division was not possible to do at 
this time because of the memory limitations of the HP-9825. 
Tais must be done pixel by pixel since division by the LSF 


1S not equivalent. The memory limitation was also the 


EZ 





reason why the HP-9825 used the LSF instead of pixel by 
exe! SUbtraction for this process. 

All the peripheral equipment to the HP-9825 was 
controlled through an IEEE standard 8 bit interface bus. 
This allowed control of all aspects of data acquisition and 
processing to be modified by software. The video recorders 
were the exception to this and provided the only real manual 
manipulations required after the program had begun. Each 
frame of video needed by the program had to be selected when 
requested by the program because of the lack of interface 
@entrol. 

The linearity of the recording system was checked and 
found to be almost distortion free. A signal from an image 
was passed directly to an oscilloscope from the vidicon and 
compared side by side to a Signal from the same image which 
had been recorded in the tape recorder, transferred to the 
video disc and then passed through the DS-30 circuitry to a 
digital analog converter for viewing. Crittenden et. al. 
(Ref. 7] present photographs of the two images for 
Semparison and further discusSion. 

Variable usage became a problem as the program grew, 
which was due to the limitations of the HP-9825. The 
HP-9825 has available twenty-six variables (A-Z), twenty-six 
Meeraoles for arrays (A=Z), a subscripted variable (r) with 


aS many elements as the memory has space for plus a 


is 





subprogram variable (p) for use in passing parameters 
through "called" subprograms. One of the usage problems was 
using identifiable variables in passing from main programs 
to subprograms. A good number of the subprograms originally 
were designed for other programs and the variables often 
conflicted. The r variable is nice to use here because of 
mess Oracvtically unlimited numbers. A more serious problen, 
though, was the inability of HPL, Hewlett Packard's version 
of the Basic language, to allow an entire array to be used 
aS a parameter to pass to a subprogram. This at first made 
necessary extra internal storage for additional arrays to be 
used only within a certain subprogram, with the variable 
name being changed through a lengthy process upon entering 
or leaving that subprogram. A much faster and easier way 
Was uSing separate storage files on the dise for each major 
array manipulation. When entry into a subprogram was needed 
the array could be stored in the file assigned to that 
subprogram, which could then be read out and labelled as any 
array needed for each subprogram. 

The programming to use the DS=-30 also was primarily a 
language problem. The DS~30 was very sensitive to the 
format in which it received data or commands. All leading 
zeroes and carriage return/line feed commands had to be 
Suppressed. The address of command locations within the 


DS-30 also had to be encoded in hexadecimal nibbles (one- 


14 





half byte) and sent in reverse order to load the registers 
properly. Commands to the DS-30 were sent as ASCII 
character strings. Data was output from the DS=-30 to the 
MEE woe> Via a fast read buffer which allowed a very fast 
data transfer rate. In use, a block of 256 bytes at a time 
were transferred representing one horizontal TV line of the 
quartered memory space. This was done in a loop which went 
back and took 256 lines, completely transferring the image. 
It was decided that output would be on the HP-9825 
plotter so that direct visual comparisons could be made 
between measured and predicted intensity LSF's. This was 
much more dramatic and easier to see than tabular output 
Which then would have to be compared on a point by point 


besos for a correlation. 
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ioe CONCEUSTON 


This program can be useful in solving the problem that 
spurred this study. It will be able to show how much wander 
can be expected out of a laser designator beam as it travels 
through the atmosphere from the laser to the target. Based 
on the video recording of the actual image on a target it 
Will show comparatively the amount of wander actually 
Present in the laser. Plans are to unfold the results 
mma ycically, but at present, a visual comparison of the 
measured and predicted patterns is made to determine if 
there is any excess wander or "jitter" attributable to the 
Mmayser platiorm itself. This determination will be helpful 
in setting specifications for the stability needed in the 
system to remove as much excessive spot movement as 
possible. It will also allow testing to determine if those 
Specifications are being met. 

This program can also be used to verify atmospheric 
problems when both the laser and target spot are fixed and 
Known. Essentially, any of the laser to target elements can 
be determined by knowing the others and calculating the 
unknown uSing adaptations of the program. The beginning 
part of this program can be used to digitize and store any 


video image and produce an LSF of the image. This may be 
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useful in recording and testing the output of several 
different lasers and optical aperatures for comparative 
ai@alysSis. 

This computer program was designed in the initial stages 
of a continuing project studying the described problem. If 
desired, changes can be made to the program to adapt to 
changing field conditions and analysis requirements. [In 
particular, the program now requires a greater understanding 
of its internal operations and operator interaction than is 
necessarily needed. More complete cueing prompts and data 
entry parameters may be worthwhile changes. Adaptability 
to other hardware may provide more program control, especi- 
ally if control of the video dise or tape can be handled 
wameeugh an interface bus instead of manually. Also, 
improvements in ocutput design can be foreseen to expand its 
usefulness and adaptability. 

Finally, the speed of program operation could be 
increased through use of a faster computer with a greater 
memory capacity. The HP-~9825 is comparatively slow and has 
insufficient memory. The use of a faster system such as the 


HP-1000 is recommended for actual measurements. 
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APPENDIX A 


USING THE COMPUTER PROGRAM 


The analysis program requires considerable preliminary 
set up and operator interaction. Before the program is 
Started there should be several items already recorded on 
videotape. The first item is a TV camera shot of the laser 
output. This is normally taken with a long focal length 
lens and filters to reduce the intensity without losing the 
beam pattern. A part of this recording should include a 
Known diffraction grating through which the beam passes in 
order to measure the scale of the data for each system of 
lenses used. The diffraction pattern should extend in the 
vertical as seen on a TV screen for proper output on the 
plotter. Also needed to be recorded is, of course, a 
sequence of laser spot images on the target. Part of this 
recording should have the target only as a background 
reference for later calculations. 

The next step that needs to be taken is to compute and 


aN then 


record the square atmospheric turbulence index Ge, 
record the following system data for later program entry; 
Wavelength of the laser, diameter of the optics objective 
lens, ratio of the obscuration to the diameter of the 


objective lens, distance to the target, and the extinction 


coefficient. 
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The first cue the program asks is for the input/output 
parameters. Table 1 lists the parameters and the different 
function each one is used for. 

If the laser source pattern is not to be computed the 
next cue requests that the laser pattern be displayed from 
the video disc for transfer to HP-9825 memory. A cue is 
also presented asking for the horizontal width of the images 
to be taken from the DS-30. The first value should be the 
left most pixel desired, the second being the furthest right 
pixel to be used. The first value must be less than the 
second, with the first no less than 1 and the second no more 
than 256. This cue is also presented below before either 
the scale image or the target image is taken. 

If the scale of the source data is to be determined the 
next cue requests that a frame be displayed on the video 
Seeoee or tne laser output with the diffraction grating in 
place. When this has been done continue the program and it 
Will plot out an LSF of that image whereby the scale of the 
data can be determined in microradians per point. 

The next cue alerts the operator to ensure that a frame 
of background video is displayed from the video disc for use 
in the difference process. 

The program will then ask for the number of total frames 
that the user wishes to average for the measured target spot 


LSF. The value should be in the range from 1 to 600. The 


Ie) 





latter value being the greatest number of frames the video 
disc can hold at one time. The program will then proceed to 
record and process images from the video dise and DS-30. 

The program will provide a cue when it wants the next image 
to be sequenced on the video disc. 

The program will alert the operator when it is about to 
produce a plot. If the HP-9862 plotter is not ready simply 
Stop the program until a piece of plotting paper is in 
place. 

The program will request the system data after the 
measured image spot data has been plotted. The data format 


mommOou Critical but should be in the correct units as 


moll OwWSs : 
Wavelength Meters 
Diameter of objective lens Meters 
Obscuration/objective None 
Beale of data Microrads per point 
Range to target Meters 
a= Meters7</ > 
Extinction coefficient Inverse meters 


If the source pattern is to be computed the program 
requests peak amplitude and standard deviation of the 


mreensity distribution. 
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The program will then continue computations without 
muicnner Inputs Until 1t has run through with the initial 
parameters. The I/O parameters will determine at which 
pmeps OUtpuUtT is plotted. After one run through, the program 
Will request a new oe range, and coefficient. If the 
input range is positive the program will again run through 
from the atmospheric calculations with the new range. If 
the range is negative, the branch is to the beginning of the 


prediction phase where all new system inputs are requested. 


Hietne range input is zero, the program ends. 
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APPENDIX B 


PROGRAM VERIFICATION AND FLOWCHARTS 


I. MEASUREMENT PHASE 

Section I (steps 0 = 50) consisted of program setup 
information and the DS-30 function command programming. 
This section was easily verified by observing that all 
initializing functions were accomplished and the response 
of the DS-30 could be visually noted. 

Section II (steps 51 - 86) consisted of transferring 
the contents of the four quarters of the DS=-30'S memory to 
the HP-9825 and producing an LSF of each image. For 
testing, an image consisting of all black (8 bit binary 
ones) or all white (8 bit binary zeroes) was placed in each 
of the four quarters of memory, each quadrant was then read 
out to the HP-9825. The values received by the HP-9825 
corresponded to the contents of the DS-30 memory quadrant 
and did show proper changes when the subject quadrant 
memory was inverted from black to white. The LSF was 
produced by averaging the pixels in each horizontal line of 
video, each line producing a point value as the image was 
scanned vertically from top to bottom. Examples of LSF's 
produced for a "black" and a "white" image are in Figures 2 


and 3. 
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section III (steps 87 - 118) subtracted a reference 
background from an image to yield a difference image, then 
shifted the resulting LSF so that the center of area under 
each image was centered on zero. In this way the wander 
from each succeeding image could be determined, removed, 
and the total frames used could then be averaged to show 
the short term atmospheric action on the laser spot image. 
Proper operation of this section can be seen in Figure 4 
where two images of slightly different position and 
intensity were averaged to yield an LSF halfway between the 
two image functions as expected. Figure 5 shows a differ- 
enced LSF where a small laser spot was added to a reference 
background. The result is almost entirely the laser spot 


With a small amount of "video noise" caused by the camera 


and associated electronics. 
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Peewee REDICTION PHASE 

Section IV (steps 119 - 139) enters the parameters of 
the atmosphere and optics needed to model and predict the 
laser spot pattern on the target. This section also yields 
the calculated Gaussian laser intensity output if that data 
1S not measured in the next section. The output of this 
section is shown in Figure 7. 

Section V (steps 140 - 172) consists mainly of 
branching, data reading, and subroutine action. These subd- 
routines, branching and data input will be verified with 
proper operation of the program as a whole. 

Spoewiomn Vi (steps 173 - 183) calculates the diffraction 
memaved Optics effect on the laser output and produces a 
meatal point spread function of the intensity. Tne output 
menrea point source should be an airy function as results 
Sa@ow in Figure 8. 

Section VII (steps 184 - 204) again consists mainly of 
moorosrams and also calculates the product of the calcu- 
Baved optics and laser output transforms. The output is in 
Figure 9. 

Section VIII (steps 205 - 220) computes the MTF of the 
atmosphere based on the inputs from Section IV. The output 
Homan Figure 10. 

Section IX (steps 221 to the end of the main program) 


iS mainly subprogram action in computing the final LSF 
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Beedicted on tne target uSing the source, optics and 
atmospheric MTF's. It also provides a link to re-enter the 
program for changed parameters. The output of this section 


is listed under the subprograms used. 
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Peis SoUBPROGRAMS 

section X (subprogram ABEL) produces a point spread 
munetion (PSF) from the LSF computed for the target laser 
Boot. tne input and resulting output is shown in Figures 
mand 13. 

Section XI (subprogram LSF) enables the program to 
compute a line spread function for later conversion to a 
faceewhnen the input is calculated by Sections IV and VI as 
point spread functions or when the data is input in 
Beecion Vin the point spread function format. Sample 
mieue and Output is in Figures 14 and 15. 

Section XII (subprogram B) calculates the fraction of 
the power inside a circle of a chosen radius about the 
Pencter of the laSer spot. This is helpful in determining 
if the spot is intense enough on the target. The input and 
Smiepuc are Figures 16 and 17. 

Section XIII (subprogram FXFORM) is a discrete Fourier 
meanistorm program that functions both as direct and inverse 
depending on the flag set in the main programs. Several 
Fourier pairs were checked to ensure correct calculations, 
oeeegures 18 through 21. 

Section XIV (subprogram SOURCE) digitizes an image of 
@aeelaser source for inclusion in the prediction 
Sewculations. it iS a storage of DS-30 video from the 


Vidicon camera input. 
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Section XV (subprogram SCALE) digitizes an image with a 
Seaudard diffraction grating for calculating the scale of 
the data. Figure 22 is an example of this image. 

Section XVI (subprogram OUTPUT) yields in useable form 
all the data required for analysis. This subprogram was 
Meedq for all the output in the plot format. Functions 
"MIN" and "MAX" are used by the “OUTPUT" subprogram. 

The program was run as a whole numerous times with 
different I/O parameters and data in order to exercise the 
branching combinations and program steps. Table 1 lists 
the I/O parameters available and their uses. This 
operation was simply to ensure the program ran as was 


intended. 
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ARRAY VARIABLE 


0¢1) 


Oe) 


0(¢3) 


O(4) 


0(5) 


TAS et 


Peo OUR UTS CARAMETERS 


VALUES 


0 


20 


0 
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BLS) a 


Same as 1 except program after one 
time through ends, not allowing.for 
computations with new range, o : 
ext. coeff. 


SOUReCe Dac vern INnDUL 1S point 
spread function. 


source pattern input is line spread 
fUMerZon. 


Source pattern input is modulation 
teansrer tuneccion, of laser only. 


Source pattern input is modulation 
Beast ehmrunietion of source and 
Opules. 

Source pattern is to be calculated 
with amplitude and sigma of inten- 
Sie yaarsre a ouulon as factors. 
Somehesmattieraer Ton anput aS PSt . 
Oortesedtfurder10n imout as LSF . 


Speres Glitiracvi1on Dattern to be 
calculated as PSF. 


No plot of source line spread 
PUneUNMOneLs CequiLred . 


Plots source line spread function 
Cavenelapleroniy if O001)=05 1 oft). 


No plot of source MTF required. 


Plots source MTF (available only if 
OCmaces new equal 3 or 4). 


NOmOoUGEmol DGEOducl of Source and 
Opures Mir required. 


DiGtsmDnOGuect mel Source and optics 
MTF's (available only if 0(1) does 
not equal 4). 





3) 





TABLE 1 (CONT'D) 


ARRAY VARIABLE VALUE USE 
Co) 0 NewplOumommopcics diffraction PSF 
required. 
>0 PlG@mewopenes ditfraction PSF 


(available only if 0(1) does not 
equal 4 and 0(2)=2). 


0(7) 0 NOmpEGemenmoptitoes ESE pequired. 


>0 Plots optics LSF (available only if 
O0(1) does not equal 4 and 0(2) does 
not equal 1). 








0(8) 0 No plot of optics MTF required. 
>0 Plows Opcres MIF Cavallable only if 
0(1) does not equal 4 
oe) 0 No plot of atmosphere MTF required. 
228) Plots atmosphere MTF. 
01710) 0 No plot of system MTF. 
>0 Plots total system MTF. 
i011) 0 No plot of LSF of predicted target 
Spo. 
>0 ElOusmlor sor  targeu Spot. 
iGl2.) 0 No plot of target PSF required. 
>0 Plots target PSF. 
0(13) 0 No plot of the fraction of power 


inside circle required. 


>0 Plots fraction of power inside a 
Ciipememot —Tadius tf. 


0(14) 0 No plot of calculated source PSF 
required. 
any 
other Plots calculated source PSF. 
value 
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abies | (CONT'D) 


ARRAY VARIABLE VALUES 
0(15) 0 
>0 


USS 
Plots measured LSF of target spot. 
Usedmas indicators of output 


branching only, entry value not 
SiQmiLF icant. 


a GP gh Sg A gp Pe fa St 


B16) 1 


any 
other 
value 


Mees msource With diffraction 
grating from DS-30 for scale compu- 
bacromsetrom OoutpUuL. No output 
DLOt ced. 


Same as 1 except plots the above. 


Does not do above input 
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AP PENDIA C 


PROGRAM LISTING 


Perens 25) ,c515),5217]) ,f3(272),0$(5),PS1(5] ,Qa{5],x5 (40) ,R 
peee) ,1.{512) 

i cd im etd Ne Wes), S140), 20140) 

meetiles O,T,£,8,G,V,W,C,1,F 

Seeermt it,Zz,c 

4: time 150d 

bree ng 3;cftg 2;ctg l 

meeread 7,l:sprt 7,R[*|];rread 7,1 

feerout “DATA",IS, 3 

Bemeecead 4,l;sprt 4,R(*) 

9: beep 

moe INPUT I/O CONTROL O[I]j" 

iieeerOor I=l1 to l6é;z;ent O[I]; next I 

fete Oll}=2;gso “SOURCE” 

ieee Cir /04 

mee tt OflojJ=l or Of[16)=2;9sb "SCALE" 

15: beep;dsp “ENSURE VIDEODISC OWN THEN PRESS CONTINUE" 
ioe st 

i7: Broo vent "# OF FRAMES TO AVERAGE", H 

ieeeeent WINDOW WIDTH ,FIRS? PIXEL",r9, "LAST PIAEL",clu 
19: H+tl-dH 

20: 07r3 

migeecOr O=]1 to H by 4 

PeeeeeOr L=2 to 9 

23: chnar(47)&cnar (55) &cnar. (64) &char (80) &char.(94)+CS 
Meeeeeit I=2;"930 VU20"+dDS 

Poetic L=3;"020"sDS$ 

Pomeett [=4;"930 L20"*DS 

Meee lt [=5;:"020"+DS 

Moeeeeit 1=6;"930 220"s0S 

@oeeeit [=7-"020"-+DS 

ewes lf I=8:"930 320"+DS 

Bees if I=9;"020"4nDS 

Beeeesc len(DS)<4;gto 38 

Peeecilar (1) &* 2"&char (2) &CS&cnar (15) &DS[5,7]*ES 
fee Cli 7 

Pee walt 100 

Som wet 704,BS 

oye cli 7 

Somerenar (46) +CS$. [1,1] 

Peemecnar{])&"2 “schar(2) &CS&char (15) &DS(1,3]*+AS 
202 cli 7 

dl: wait 100 

42: wrt 704,AS 

eee Cli 7 

44: if [#3 or O>1;gto 47 


oye, 





dsp “BACKGROUND RECORDED" 

walt 5UU0;dsp "CONTINUE WHEN NXT FR de READY"; St 
te DSE"020" 0 or I=3 or I=9 ana Ox 57g to 

asp “CONT. WHEN NXT FRAME ScLECTED  ;stp 

next I 

Mme OrLha4;qtO 53 


: for X=l1 to 4 


gto 54 

hor A=l €O H=-O+1 

ue I Loa Cg ketot DP oe har: (96)+PS 
Pex g i; gto 5/ 

l+J;0+Y;gto 62 

Memxt es GEO 59 

Bug; 22070; 0>¥;gto o2 

Pee ato; gto ol 

ed 2 55eY 7; GEO O2 
Oe 32 55+Y 


: for L=l to 256 


asp L 

O>D;3""s0S; Y+1>Y 

Y*512+J+S 

for I=2 to 6 

Poti o/ 16) +T 

T>U 

S-T*16+1I* 16-0 

Uses 

OS&cCnar (.T) +08 

next I 

char (1) &"="&char (34) &char. (2) &0S&char. (3) &P S&char.(15)+5S 
mit 1,2, Cc 

oii 7 

wrt /04.1,BS$ 

put “pATA" 

met /04, “DATA” ,256 

mast DATA™)+E;i1£ E=-L;jmp 0 

for I=r9 to rl0Q 

num(IS.{I,1I ]) +2 

D+Z2D 

nexc If 

D/ (r10~-r9) >R[L] 

next L 

MeL lg3=0;gto 91 

"TAKE DIFFERENCE OF IMAGE AND BACKGROUND": 
Mmeeaqd 3,1>sread 3,1(*}];£or I=l1 to 256 
abs(R[I]-I[1I])*+R{[I];next I 

Geo 93 

meeagd) 3,1;sprt 3,R[ *] 


eeageto 113 


UeG :07F 
per L=30 to 226 
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95: R[I]+GrGrnext I 

poet Oor 1-30 to 226;R[I1]+F>F 

Bye lf FO=G/2;gto 99 

98: next I 

99: [+A 

mee: if O=l and X=2;lerl 

lull: rl-A+r2 

Me fOr J=l to 250tr2 

103: R{J-r2] +R[J] 

mee: 1f£ J<256tr2:next J 

ieee fOr K=25o+r2 to 256 

100: OFR{[K] 

igy: next « 

Mee: Sport /,r2 

109: rread 4,l;sread 4,1[*] 

mee: £Or I=l1 to 256 

lll: I{I]+R{I]+R[I] ;next I 

Mmm: rread 4,il:sprt 4,R[*] 

Meee next A;next O 

114: "TAKE AVERAGE OF ALL FRAMES": 

Meee EOr t=) to 256; R{[1]/(H-1) +R[I];R{[1]+R[513-1I]snext I 
eos creaaq 4,l;sprt 4,xR[*] 

ieee: 1t O[15)]<l;gsb “OUTPUT” 

118: beep 

Meee ent "RANGS",R,"SQCN",O,"EXT CCEF",£, "WAVLNGTH",w,"OBJ LE 
NS" ,O 

mmpeent OfG/0BJ , 58, "SCALE",rll;:sprt 6,R8,0,£,w,0,8,rll 
121: "SYSTEM DATA-W=wWAVLNTH IN METERS O=DIA OF OBJ LENS IN ME 
wR Ss" s 

122: “B=DIA OF OBSUR/DIA OF OBJ LENS rll=SCALE OF DATA IN AIC 
RORADS/PT": 

fess) R=DISTANCE TO TGT IN METERS": 

124: "Q=SQUARE OF ATMOS TURB INDEX IN METERS” -2/3": 
Mes | })6vB=—EXTINCTION COEFF (1/METERS) ": 

mee: tf O[lJ#5;gto 14u 

127: beep 

Wee. ent “AMPLITUDE” ,A, "SIGHA",C 

fee £Or I=l to 250 

130: rll*(I-1)+F 

ieee FF 2/(2*C" 2) +G 

Meee LE G>i3;gto 134 

ies; A~exo(-G)erI./1];1+K; gto 135 

mea; O+I[T1} 

iy: next I 

moe crread 5,l;sprt 5,1[*]:rread 5,l1;sread 5,R[*] 

meee L>O(15);if Of14]#0;9gsb “OUTPUT" 

ioe: gto 149 

meg: “SOURCE PATTERN DATA INPUT": 

140: rread 2,l;sread 2,R[*] 


Meee £Or I=l to 256 
142: I-+K 


a7 





meee Lf Ri 1}<—500;gto 146 
ive: next I 

mes: gto 149 

Mae: tor J=I to 256 

ag? UR J] 

148: next J 

149: O{1lJ+I 

P5003; “BRANCH ACCORDING TO TYPE OF SOURCE PATTERN DATA": 
omceomlaa tito ie gtO 155-1 I>2:;gto lol;irt I>3;gto 206;1f 
Mee oto 152 

152: "CONVERTS SOURCe POINT SPREAD rUNCTION TO LING SPREAD FC 
Ne 

leer: 9GSb “LSF" 

Meee 2>O(15);1£ O13) >0:;gsb “OUTPUT” 

155: “CALCULATION OF FOURIER TRANSFORM OF LINE SPREAD FCN OF 
SeuKCE”™ : 

foo: fOr I=l1 ts 256 

ieee Ri L|+kK{[S5S13—-1] snext I 

ies 9Sb “PF XFORM" 

fees Cread 5,1l;sprt 5,R[*] 

Meee SeO(15);1£ Of 4] >0;gsbp “OUTPUT” 

meme: if Of 2) =2;gto 173 

ies OPTICS DIFFRACTION INPUI Ir NOT CALCULATED": 

lo3: rread 8,1;Sread 38,R[*] 

iver. £LOr I=l to 256 

lo5: I-+«K 

Mages LE R[I)]<-500;qto 169 

Moy: next I 

Mee. gto 182 

169: £o0r J=I to 256 

fee USR(J)snext J 

mee gto 182 

fee CALCULATION OF DIFFRACTION LIMIT POINT SPREAD FCN": 
173: B° 24D; 1-D-H; 3.14159e-6* rll*O/w-Z 

174: for I=l1 to 256 

Meee 4~*(1-1)*+¥;1f Y>30;gto 178 

ees ( ALRY (Y¥)=D* AIRY’ (Y*B))°2/H° 2+R[I] 

ie: I=Ksgto 179 

mo: U>R([I] 

ioe next I 

ioe Creag 8,l;sprt 8,R{*!] 

ie S2°O(15);71F Of 6)>0s;asb "OUTPUT" 

Mere if Of 2)=l;gto 132 

183: "CONVERTS OPTICS POINT SPREAD FCN TO LINE SPREAD FCN": 
184: rread 8,l;sread 3,R[*] 

Mee: OSD “LSF"” 

Meee rread &,l1;sprt ¢,r[*] 

Mowe SeO(15)5;i1f£ Of 7) >0;gsp "OUTPUT" 

188: "TAKES FOURIER TRANSFORM TO GET OPTICS MODULATION TRANS F 
Peach” : 
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loys: 


193: 


198: 
o> : 
Zw) : 
201: 
eu2z: 
Zi0> : 
204: 
205: 
ao : 
euy: 
Z06 : 
Zero : 
Zar : 
Zeel : 
Zan : 
25, 
214: 
Zr: 
Zao ¢ 
way: 
28: 
Eas : 
Z20 : 
Ze: 
el: 
“e383: 
24 : 
Oe: 
ZO : 
Za]: 
ZZ 6 3 
ae : 


Ze: 
Z3ol: 
252: 


mecmGd earl; Sreag 6,h| * | 

POneel=1 "EO 250 

R{I]+R[513-I] ;next I 

SSbre bse ORM 

EOLte i=l tO 256 

Pon (ei gto 196 

R[IJ/R{LJ+R{[I] ;next I 

BEOdG sop, sprt G,R|*} 

SeOieojein Ole} >0;gsb “OUTPUT” 

"CALCULATION Or THE PRODUCT OF TWO FOURIER TRANSFORMS": 
rread 8,l;:sread 6,R[({*];rread 5,1l1;sreaaq 5,I1[*] 
zEOr [=] to 256 

mies LL ri>R(Ils-next IL 

magecaci0, 1; sort LO,R{* |} 

POS j-2b is} >0:-gsb “OUTPUT” 

exp(-E*R) »X 

"ATMOS AODULATION TRANSFER FCN CALC-SHCRT TERM"; 
W*R*21.6eA;3;W (-. 3333333) 3B; OV; 0+L; X*R[1] 

for I=2 to 250 

MeQuz4e ir ise -U>gtO 215;1£ L=U3sgto 209 


stg l4 

1£f Go13;1+L 

1£ G-V<0;-1+L 

ee exo (—Gioki Pf) piekK;ctrg 14;gqto Zid 

Weil, gto 216 

R({K] *cos.((I-K)*3.14159/K)°~2+R{T] 

Sire 7,18 I-3*K/2<0;gto 212 

IeKk ;1>L 

next I 

mecad 9,l7-Sprt 9,R{(*] 

peeieloje tt O[9) >O0s;gsb “OUTPUT 

"TRANSFER FCN OF SOURCE*OPTICS*ATMCSPHERE": 
rread 10,l;sread 10,R[*];rread 9,l;sread 9,i[*] 
Por tL=] to 256 

iio | Lj) >ii 1} 

ie | >R{ t } 

R{IJ+R[513-[1] 

next [ 

Peewee tr eo (lO) >0;qsb “OUT PUT" 

"INVERSE FOURIER TRANSFORM GIVES TARGET LINE SPREAD FCN" 


Saegel;9so “EXFORM"’;1+K;ctg 1 
MerOotmioj7it O(lij>O-qsb “OUTPUT 
"CONVERTS LINE SPREAD FCN TO POINT SPREAD FCN BY ABEL TER 


ANSFORM" ; 


mes: 
es4: 
Be 5: 


gsb “ABEL" 
Peres rie Oli2)/>0;gsb “OUTPUT” 
“CALCULATES THE FRACTION OF THE POWER PNSiDE, A CIRCEE Or 


map rUS RR": 


Sas, 





Zee: gsb “B" 

Mewes t2Z-O)1S}sif O[13]>6;qsb “OUTPUT” 

238: “READ IN NEW RANGE AND Q If RANGE NEG. READ IN ALL NEW D 
AA” : 

239: "STARTING wITH I/O CONTROL CARD,IF RANGE=G6 PROGRAM STOPS 


240: “IF RANGE IS POS. THE PROGRAM CARRIES OUT CALC’S WITH .TH 
bE NEW": 

241: “RANGE AND QO": 

Mee: if O[1J=O0;gto 245 

Peewent “RANGER ,"SOCN",)),"EXT COEF" ,E 

244: gto 200;1f R<=0; gto 119;1if R=O0;gto 245 

feat 6end 

meee ABEL :R([LJ*N;1.4R[LJ—-1L.8R[Z}]+.4R[3]eR[1L]}];:ror I=2 to K 
247: R[Ijem 

248: .4N+.2M-.6R[(I+1L]+R[I] seoleNsnext I 

249: for I=l to K 

Zoo: RLL)/(2*y¥( (1+.1)° 2-I1)) +R[TI] 

Zeer £Or J=I+l to «K 

252: R[IJ+R[J]/yV¥( (J+.1)° 2-11) >R[TI] 

fee next J;R(I]/m>R(I J] :dso I;next I 

254: ret 

meee «6 LSE“ :sfor I=l to K 

feo: i+Jsdsp I;R{[I] +z 

257: y (I*I+T*J) oy 

Meee 2~ ((l-frco(Y))*R[int(Y)jJ+frcol(yY)*R[int(¥)+1])}+2+Z 
gees UtleJ;if Y<K;jmp -2 

260: Z*R[I]:nexct I:ret 

Bomes BY :.250TR[LJ]>R[1] 

mee: £Or [=2 to K 

263: 2nTIR[I]+R[ IL-1] +R[I] 

264: next I 

205: 1f R[K]=U0;ret 

mmo LOr I[=l1 to K-R[1L.jJ/R[ KJ *R[I]snext I;ret 

fom 6 ALRY":if pl<0;beep;dasp “ERROR pl<O";stp 
Moore if pl=O0;ler4;ret r4 

eee Or S71f pl>l5;jmp 2 

ZU : Z0+l0pl—pl (2/3)+r6;jmp 2 

271: 90+pl/2-+r6 

Meees TE D1<5;6+pl+rl2;jmp 2 

oye 1. 4*clto0/pl->rl2 

274: max(int(rl2),int(3+pl/4) )+rl2 

275: for M=rl2 to r6 by 3;1le-28+r8;0+rl3+rl4 
Pyee- SLOG 1O0;1f M/2=int(M/2) ;cfg 10 

Z77: Lor J=l to M—-2;2* (M-J) *r8/pl—-ri3erl15;r8-rl3 
278: rl5er8;if M-J-2=U;rl5-r4 

eee cme 10; r1i4+2*r8*flglOe*rl4;next J 

Zee: 2*r8/pl—-rl3rl5 

Pees ElLAt+rl5er1l4;r4/rl4+r4 

232: 1f abs(r4-r5)-abs(r4*le-6) <=0;2*r4/pler4;ret r4 
283: r4+rS;:next M 
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234: 
Za): 
Zoo : 
Zo: 
Mmreres=0) to 2°(M-l)-lLscll BI (J,2,N-1) 


Se ° 


beep;dsp “Accuracy not obtained";wait luOO;ret r4 
"FE XFORM" :rad;92N;ina I 

1/2~ (N-1) +T 

for M=l to N; 2° (N-i) >rl6 


cos (PT) +C3sin(PT) (1l-~2f1gl) +P 


meoQretl=2r Loutl co 2ZrloJ+r1lo 


R[{[I])er0;R[I+rl6Jerl 

I{(IJer2;I[I+trl6j+r3 
rOQ+rl*C+r3*P+R[ I) ;r2+#r3*C-rl*P-I [IJ 
rO-rl*C-r3*P+R[I+r16] 

r2-r3 *C+tr1l*P+I[(I+rl6] 

mexc IT e-nexte Jsnext M 

momel=-Oubo 2 N-)lscli BI (I,J,N);if I-J>0;jmp 8 
be =I jmp 5 


: R[I+l]/y¥(2°N)>P 


I.{I+l] /y¥(2°N) +Z 
R(J+t1ljJ+R[(I+1l}) ;I([JrtlJ+>+LlLl[irl] 


>: P+R[J+1] 3221 [(J+1] 


R[I+l]/y¥(2°N) eR[I+1] 

I({I+tl]/¥(2°N) +I[1I+1] 

next I;deg;ret 

wee SUD 27 pls+p4 

por 2Z=] to D3 

p4/2+p4;2p2+p2 

Moire (p44) 470; p2+1+p2 

int(p4) +p4;next Z;ret 

BO@URCE”™ ssig 2 

beep;dsp “DISPLAY SOURCE IMAGE ON VID DISC";wait 3000 
et) “SCALE ~ 

mecagd 2,l;sprt 2,R( *} 

Seg 2 ret 

poeAltie :1£ flg2;qto 318 

beep;dsp “DISPLAY SCALE IMAGE ON VID DISC";wait 3000 
Poop tHEN PRESS CONTINUE";stp 

Pee WiibOW wilorH, FIRST PIXEL”,r9,"*LAST PIXEL",r10 


metor 1=2 to 3 


char(47)&cnar (55) &char (64) &char (8Q)&char (96)+CS 
poeet=2;°930 O20"+DS 

1f [=3;"020"+DS 

Pete n(DS)<4; gmp 5 

Seuligeilye czceecharn (2) &CSéchar (15) &D$[5,7]+B$ 
ert / 

walt 100 

wrt 704,BS 

Cnar.(46)+CS([1,1] 

Stamnes 2  scehar (2) &CSé&char (15) &DS[1,.3]+AS 
eli. 7 

walt 100 

wrt 704,A$S 
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S5u : 
So: 
Bir : 
oS: 
354: 
soo: 


So 2 
S56 ° 
B9 : 
300: 
3o0l: 
Bo : 
Bioo: 
304: 
36» >: 
366: 
BO : 
S06 : 
369: 
370: 
Se 
Bye: 
3 is 
274: 
m5: 
B70 : 
ey: 
Bo: 
379: 
SoU : 
Bol: 
Bio : 
333: 


next I 

AP Se a ugar. »  echar (96)>+PS 
l+J;0+Y 

BoE L=1 to 250 

dsp L 

DeDs ""sSOSsYt1>Y 

Y*512+J+S 

fon. [L=Z. CO.6 

int(S/16) +T:T+U;S—-T*16+1*16-T 
Us*S;CS&cnar (.T) +O$S 


> hext I 


char (1L)&"="&char (34) &char (2) &0S&char (3) &PS&char (15) +BS$S 
eli 7 


mewet 7U4.1,B5 


burt "DATA" 

mee 704. DATA", 256 

rds ("DATA") +E;1if£ E=-l;jmp 0 
for I=r9 to rl10 
num(IS[I,I])+Z;D+2+D 
next I 

D/ (c10-r9)+R[L] 

mexc. L 

Per elg2;gqto 358 

MEO |lo|)=2-9gsb “OUTPUT” 
ret 

Oupe wT ” 3 

beep 


asp "STOP IF PLOTTER NOT READY THEN CONTINUE";wait 5000 
stg 14; ‘MIN’ (256) / “MAX°(256)+rl7;cfg 14 

Pe Coeiesecl —,75,1.05,—-30,26043plt 1,0,1; jmp 2 
Soe eetoyt.05,-30,264-plt 1,0,1 

momel=!0) to 0 by —-l 

primey LO,0,2;7plt 1710 ,256/150,2 

pare wl/ 10,0, 2:;next I 

Eon L=0 to 16 

emmy eltol, 2,01 1/7150,161 ,2 

mie 0,161 ,2 

next I 

Pemrosiz .9,1,1.4 

mem se Ome l=167to 0 by -~L-plt -.1,161-1,1 

1lbl 16(16-1) ;next I 

Pacer LO l—-0 to 1Osplit 1/10-.06,-5,1 

Webel /lue-next I 

1£ Of15J=10 or O[15]=0;"TARGET LINE NUMBER"?YS 
MeO 5)=1 or Of15j=11;" POSITION VALUE"+YS 

1f O[15jJ=10;"PREDICTED NORMALIZED LINE VALUE" +xXS 
if O[15]=0;"MEASURED NORMALIZED LINE VALUE" +xX$ 
Pate oie -eSlzZ f.2,1,1.4,90;1lbl YS 
MipcwUm=iS,1°eSiz 1:2,1,1.4;1bl «S$ 

SoZ ply « f > OU 
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384: “MAX* (2560) >r7 

Boos) Lic 7=0; Ler? 

Seo: £Or J=l to 250 

Meee Olt Ri J) /r7,257-J;next J 

388: pen;ret 

bog: 'MIN“ sR[1L]>p2 

Sees £Or L[=2 to plsift R[LT]>p2:jmp 2 
Bee R{[I)|+p2 

mee: next I;ret p2 

Bese: MAX" :R({1]+p2 

Bere sor I=2 to pl;1f R({I]<p2;j3mp 2 
395: R{[Ij+p2 

396: next I;ret p2 

eel io4 
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